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Electrical transport and specific heat properties of Ndi-^Pb^MnOa single crystals for 0.15 < x < 
0.5 have been studied in low temperature regime. The resistivity in the ferromagnetic insulating 
(FMI) phase for x < 0.3 has an activated character. The dependence of the activation gap A on 
doping x has been determined and the critical concentration for the zero-temperature metal-insulator 
transition was determined as x c ~ 0.33. For a metallic sample with x — 0.42, a conventional electron- 
electron (e-e) scattering term oc T 2 is found in the low-temperature electrical resistivity, although the 
Kadowaki-Woods ratio is found to be much larger for this manganite than for a normal metal. For a 
metallic sample with x = 0.5, a resistivity minimum is observed for x = 0.5. The effect is attributed 
to weak localization and can be described by a negative T 1//2 weak-localization contribution to 
resistivity for a disordered three-dimensional electron system. The specific heat data have been 
fitted to contributions from free electrons (7), spin excitations (P3/2), lattice and a Schottky-like 
anomaly related to the rare-earth magnetism of the Nd ions. The value of 7 is larger than for normal 
metals, which is ascribed to magnetic ordering effects involving Nd. Also, the Schottky-like anomaly 
appears broadened and weakened suggesting inhomogeneous molecular fields at the Nd-sites. 

5-H 1 PACS numbers: 71.27.+a, 75.47.Gk, 72.10.Di, 65.40.-b 
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I. INTRODUCTION 

Colossal magnetoresistive manganites (Ru-^A^MnOa ) are known for interesting physical properties and compli- 
cated electronic phase diagrams [H, S S 0; I3| ■ In particular, their intrinsic properties in the low- temperature phases 
are far from trivial. While many of these properties were known for more than 50 years [6|], an appreciation of the 
magnitude of these effects is a more recent development Q- Among the different phases in these manganites the 
occurrence of a ferromagnetic insulating (FMI) phase poses a difficult problem because the double exchange model 
I I predicts only metallic ferromagnetism at low temperature [H, 0, [T0| • The ferromagnetic insulating behaviour is ob- 
■ served in manganites when the divalent dopant concentration is quite low. At higher dopant concentration, crossing 
CN| ] the threshold limit, a ferromagnetic metallic phase (FMM) is usually found. The low-temperature transport properties 
of manganites have been studied in some detail only in recent years [2J. Urushibara et al. described the temperature 
dependence of the electrical resistivity in Lai-^Sr^MnOs in the composition range of metallic conductivity at low 
temperature as p(T)= p(0) + AT 2 fill. Similar data were reported by Schiffer et al. [H| and Snyder et al. [Hj]. As 
pointed out by Kadowaki and Woods [lj| , a general relationship between the coefficient A for the T 2 contribution to 
the resistivity and the square electronic heat-capacity coefficient 7 is obeyed in usual metallic systems. However, in 
manganites very large ratios A/j 2 are found that are an order of magnitude larger than the Kadowaki-Woods ratio 
for metals with strong electron-electron (e-e) interactions [15j |. This casts doubts on the assumption that the usual 
electron-electron (e-e) scattering underlies the T 2 dependence of the resistivity observed in metallic manganites Q. 
In particular, manganites are well known examples for bad metallic behaviour pH, [l7|. 

Some recent reports find that the manganites have resistivity minima and the resistivity at low temperature is 
higher than Mott's maximum metallic resistivity of about 10mO cm [HI, EH H(| ■ According to the scaling theory for 
disordered electronic systems in 3D metals [2l[, the weak localization and electron-electron scattering in presence of 
strong disorder should influence the electrical conduction in manganites. This type of e-e interaction is described by 
a negative T 1 / 2 power law contribution to the resistivity. This contribution has been employed to fit the resistivity 
data of the intrinsically disordered systems like manganites in recent reports [l8l . Il9l . |20T | . Among other mechanisms 
ruling the temperature dependence of the transport properties, the contribution of electron-magnon scattering to the 
resistivity in manganites poses further problems. In the ground-state, one-magnon processes should not take place 
because of the half-metallic band-structure of the manganites with a fully polarized band of the electrons at the 
Fermi level (22J. Additionally, the appearance of a Kondo anomaly in this kind of ferromagnetic compounds with 
spin-disordered clusters has also been discussed [2fj| . This also suggests that there is a need of some modification in 
the dependence p(T)= p(0) + AT 2 to describe the low temperature metallic transport property of manganites. 

A new theoretical evolution, mainly triggered by a strong-coupling approach to include the electron-lattice coupling, 
can now describe the transition between the FMI and FMM with varying doping. Pai et al. [23| have recently proposed 
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a new effective low-energy Hamiltonian starting from two qualitatively different coexisting vibronic states at each site 
of the lattice these being labeled as I and b. Here / vibrons describe a localized Jahn- Teller (JT) polaron and the 
other b vibrons form a broad and dispersive band. Within this approach, the insulating gap in the FMI phase at 
low temperature is identified as the T — electrical gap A between the occupied I levels and the unoccupied b 
band bottom. In principle, the theory is able to make detailed predictions on the electronic structure and transport 
properties through the whole range of compositions between the FMI and FMM phase in the perovskite manganites. 

These theoretical developments and other experimental findings inspired us to explore the low temperature transport 
and specific heat properties of mixed- valent manganites. The present analysis is based on experimental data obtained 
on a series of sin gle- crystals from the Pb doped NdMnC>3 system, which is a less studied member of the perovskite 
manganite family [2J, |25j. From our own earlier studies on these crystals, we presented a tentative phase dia gram of 
Ndi-^Pb^MnOa for wide range of temperature and composition (26| and some of their physical properties [26l.l27l.l28l 
[29| . In the present investigation, we report further properties of these Ndi-^Pb^MnOs single crystals and a detailed 
analysis of their low temperature resistivity and specific heat data. 

II. EXPERIMENTS 

Single crystals of Ndi _ T Pb r MnQ3 were grown by a high temperature solution growth method using PbO/PbF 2 
as flux, as reported earlier [25| . The compositional analysis has been carried out by Energy Dispersive X-ray (EDX) 
and then substantiated by Inductively Coupled Plasma Atomic Emission Spectroscopy (ICPAES). Resistivity mea- 
surements were carried out by the standard four probe method in the range 4-300K. The electrical contacts were 
generally made on the (100) plane of the single crystal samples. The contact material was a Ag-I5 weight % In 
alloy. Magnetoresistance (MR) was measured in similar temperature ranges at 7 T magnetic field. The specific heat 
measurements were carried out in a physical property measurement system (PPMS model 6000, Quantum Design) 
from 40 K down to 2 K. 



III. RESULTS 
A. Transport 

Resistivity data for the Ndi-^Pb^MnOs single crystals were already presented in Ref. 26] focussing on the metal- 
insulator (MI) transition at intermediate temperatures around 150 K. The activated behaviour, p — po exp(_E/fcsT), 
is put into evidence by the logarithmic plot of the resistivity p as a function of inverse temperature 1/T. This is 
shown for x = 0.15, 0.3, 0.4 and 0.5 in Figure 1. The slopes of the straight sections of the resistivity in this plot 
measure the activation energy E. The crystals with x = 0.15 and 0.3 display an activated behaviour in both the 
paramagnetic and ferromagnetic state, albeit with different activation energies. The activation energy E = Ea in the 
high temperature parama gnet ic insulating phase can be related to transport through Jahn-Teller(JT) polarons, as 
discussed in the literature [301 ] . The activation energy Ea depends on x and becomes minimum for the most metallic 
sample close to x = 0.4 in our series of crystals [26| . 

For the ferromagnetic insulating samples, the same type of behaviour holds with an activation energy given by a 
gap energy, E — A, in some temperature range below Tq towards zero temperature. The gap energy A, as derived 
from activated fits, decreases with the dopant concentration x. The dependence of the gap energy A versus dopant 
concentration x is shown in Figure 2. It can be compared with theoretical predictions. The critical doping for the 
T = ferromagnetic insulator to ferromagnetic metal transition is to be determined by the vanishing of A. According 
to the theory of Pai et al. [23[ a dependence A(x) = Ejt — x 1 ! 2 Do is expected, where Ejt is the energy scale for 
the electron JT-phonon coupling, and Dq is a measure of the electronic band-width of the itinerant e g electrons of 
Mn. Fits for the doping dependence of the experimental values of A(x) to this square-root behaviour are shown in 
Figure 2. The values obtained from this fit for Ejt and Dg are around 273 (±38) and 493 (±79) meV respectively, 
when fitting was carried out in x— 0.15-0.3 range. On the other hand, the fitting in x= 0.2-0.3 range yields Ejt and 
D as around 195 (±12) and 339( ±25) meV respectively. It is observed that the system for dopant concentration x 
= 0.15 is not well fit by the theoretical square-root dependence for A(x). As seen from Figure 2 the experimentally 
determined value for A(x = 0.15) is appreciably increased. It is probable that at this low doping level further effects 
increase the apparent activation energy for the electronic transport. The extrapolation of A(x) towards a vanishing 
gap A = yields a critical doping level for the zero-temperature MI transition. The critical concentration for the zero 
temperature MI transition in this system has been calculated from the parameters obtained from the fitting excluding 
the point corresponding to x = 0.15 giving a value x c = (Ejt/ Do) 2 — 0.33(±0.02). We consider that the fit excluding 
x = 0.15 data provides the more reliable estimate of the critical concentration. 
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Now, assuming that A vanishes for this critical doping level x = x c , we can proceed to analyze the transport data 
for samples with x = 0.4, 0.42, 0.5 > x c , which display metallic character at low temperatures. Essentially, x= 0.4 
and 0.42 show similar behaviour. The resistivity at low temperature can have contributions from residual resistivity 
Po due to static defects independent on temperature and from scattering by elementary excitations. The effect of 
electron-electron scattering can be described by a term proportional to T 2 15]. Taking these two contributions into 
account the resistivity should follow the relationship as p = po + AT 2 . The data for x = 0.42 (Figure 3(a)) are 
expected to follow this conventional behaviour. However, the p(T) dependence with x— 0.5 sample shows a shallow 
but distinct minimum in the low-temperature region (see inset in Figure 1) indicating further effects. The resistivity 
minimum persists in the presence of magnetic fields and even becomes deeper, while shifting to higher temperature 
of about 110 K in a 5 Tesla field (Figure 3(b)). We assume that this anomaly is related to weak localization and 
electron-electron interaction in the correlated electron system owing to static disorder and large Coulomb interactions 
[2l| . Additionally, higher order contributions in temperature may play a role. These are two-magnon interactions 
with a term oc T 4 5 and an electron -phonon (e-p) term oc T 5 . Because it is difficult to distinguish between these two 
terms, we do not consider these two possible contributions separately, rather we restrict our fits to the two-magnon 
term only. Hence, the resistivity for the metallic low-temperature range should be described by 

p{T)=p» + AT 2 + p t T 1 l 2 + Pm T^. (1) 

For the resistivity of the x — 0.42 sample, no indication of weak localization could be observed. To verify this, we 
have plotted the resistivity data after subtracting the residual resistivity against T 1 / 2 (see inset of Figure 3(a) ). The 
data are not linear as a function of T 1 / 2 in that temperature range, whereas linearity is supposed to be a prominent 
signature in support of weak localization and a T 1 / 2 -term due to e-e interactions [U, [Hj]. Hence, resistivity data for 
x = 0.42 have been fitted according to equation (1) without the T 1 / 2 term. The result is shown in Figure 3(a). The 
fit parameters for the stable fit in the temperature range T < 50 K are listed in table 1. 

TABLE I: Summary of fitting of the resistivity data for the x — 0.42 single crystal in the low temperature range according to 
equation. (1) with error bars and standard deviations (x 2 )- 



x values 


Po 


A 


pm 


T 


x 2 




ohm- cm 


ohm- cm K -2 


ohm-cm K" 4 - 5 


K 




x =0.42 


0.2001 


7.4 xl0" b 


6.5 xl0 _il 


<50 


5x 10 _Y 




± 0.2 xlO" 3 


±4 xlO" 7 


± 2x 10" 11 







The resistivity data of the x = 0.5 crystal require inclusion of the p^T 1 ! 2 term due to weak localization. This is 
obvious from the linear part of the plots p versus T 1 / 2 in the Inset of figure 3(b) with and without applied magnetic 
field. However, the presence of the electronic contribution ^4T 2 to resistivity is uncertain. As listed in Table 2 we 
have tried various fits for the data in zero magnetic field: a fit in the low temperature range excluding the T 4 5 term 
in the temperature range T <60 K leads to an insignificantly small or even negative contribution AT 2 with large 
errors. Similar results for the AT 2 term are found for fits by including the T" 4 5 term in an extended temperature 
range T < 100 K. Therefore, this electronic term appears to be irrelevant for a valid description of the data. However, 
one can get an upper estimate of the parameter A ( 6.3 x 10 -6 ohm cm K~ 3 / 2 ) from the fit in the range T <60 K. 
By putting A = 0, we find a satisfactory fit in the same temperature range in zero field. A fit with similar quality 
is achieved for the resistivity data in a field of 5 Tesla (table 2 and figure 3(b). Here, the magnitude of po is slightly 
increased as compared to the fit to the zero-field resistivity. It is to be noted that the lowest possible temperature, 
at which we could measure the resistivity is around 5 K for x = 0.5. Hence, the value po is actually close to p§K- In 
the fits, there is a trade-off between the negative T 1 / 2 term and the po, which leads to some systematic interrelation 
between these contributions. Therefore, the shift in po is probably not related to a real physical effect due to magnetic 
field. However, the strong increase in the magnitude of p e due to the magnetic field is clear, see figure 3(b). The 
two-magnon term T 4 5 is suppressed in a magnetic field, which suggests that it is dominated by two-magnon scattering 
contributions that are subdued by the magnetic field. 

To ascertain the essential absence of the expected electronic contribution to the resistivity AT 2 , we have additionally 
used graphical methods. We plot the resistivity against T 2 (after subtraction of the residual resistivity). The data for 
x = 0.5 do not follow a straight line (inset (b), Figure 4), unlike those for x — 0.42( inset (a), figure 4). In order to 
better display the deviations from the T 2 contribution, we also have numerically differentiated the data and plotted as 
A^ 1 dp/dT 2 vs T in figure 4. Here, we have used the calculated value of A from the slopes of the straight line fits (see 
insets, figure 4). In the plot of figure 4, the data show an upward deviation from 1 at higher temperatures. This should 
be related to the expected weak two-magnon scattering ( T 4 5 ) and electron-phonon scattering (T 5 ) contributions. 
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TABLE II: Summary for fits of the resistivity data for the sample with x at noH = and 5T in different temperature ranges 
(x 2 is the standard deviation) 



x,H values 


PO 


A 




Pm 


T 


'2 

X 




ohm-cm 


ohm-cm K 


1 T7- 1/2 

ohm- cm K ' 


ohm-cm K -4 ' 5 


K 




z=0.5 


0.626 


-6.2 xlO - ' 


-1 x 10 -2 




<40 


3.4 xlO -Y 


H=0 


±0.1 x 10 -2 


±9.34 xlO" 7 


±0.3 x 10" 3 










0.6327 


6.3 xlO" 6 


-1 x 10~ 2 




<60 


2.4 x 10 




±0.1 x 10 -3 


±7.4xl0 -7 


± 0.5 xlO -3 










0.6296 


5.26 x 10" 7 


-1 x 10~ 2 


1.6xl0" 10 


<100 


2 xlO" 6 




±0.1 x 10 -3 


7.4 x 10 -7 


±0.4 xlO -3 


±6 x 10" 12 








0.6287 




-1 x 10 -2 


1.6 xlO ~ 10 


<100 


2 xlO" 6 




± 0.12x 10 -2 




±0.2 xlO -3 


±2xl0" 12 






x= 0.5 


0.6513 




-3 xlO" 2 


2.8 xlO -11 


<180 


2 xlO" 5 


H=5T 


±0.18 xlO -2 




±0.2 xlO -3 


± 2 xlO -13 







B. Specific heat 



The specific heat C vs temperature data for the series of Ndi-xPbzMnOa crystals are plotted in Figure 5. For 
manganites the basic contributions to the specific heat in the low temperature range can be described by [33| . 

c = /3 3/2 t^ 2 + jT + b 3 t 3 + b 5 t 5 . (2) 

Here, /3 3 / 2 is the coefficient of the contribution from spin wave excitations for ferromagnetic order, 7 is the coefficient 
of the electronic specific heat, B3 and B5 are coefficients of the contribution from the lattice. In the temperature range 
T < 15 K, an additional contribution from a Schottky-like effect strongly influences the behaviour of the specific heat 
for the Nd-based systems [13] ■ The best range for the fitting according to Eq. ([2]) has been found to be 20 — 40K. 
For this temperature range nuclear hyperfinc effects in the specific heat need not be considered, as they contribute 
appreciably only at much lower temperature T < 2 K. The results of the fitting by Eq. @ is shown for x — 0.4 as a 
representative plot in Figure 6. Similar fits have been achieved for x — 0.15 and 0.5 (not shown in the figure). An 
analysis of the specific heat of the a;=0.3 crystals has already been reported in Ref. |27| . 



TABLE III: Summary of fitting for the specific heat data in the temperature range 20-40K. The values of the spin-stiffness D 
derived from magnetization data via Bloch T 3 ^ 2 law have been used to fix the ^3/2 parameter for the ferromagnetic spin- wave 
contribution. The definition of the other coefficients are given in the text. \ 2 is the standard deviation. The parameters for x 
= 0.3 are taken from Ref. [13] 



Composition 




7 


B 3 


B 5 


D 


x A | 




mj mole - 


mj mole -1 


mj mole - 


mj mole -1 


meV 






K -5/2 


K -2 


K -4 


K-6 


A 2 












x 10 -3 






Ndo.ssPbo.isMnOs 


19.57 


71.82 


0.26 


-0.068 


17.3 


4.8 xlO -3 






±0.11 


±0.01 


±0.004 






Nd .7Pb .3MnO 3 


6.51 


71.73 


0.356 


-0.98 


34.9 


1 xlO -3 






±0.13 


±0.05 


±0.007 






Nd .6Pb .4MnO 3 


6 


68.64 


0.387 


-0.112 


37.8 


3.3 xlO" 3 






±0.01 


±0.004 


±0.002 






Ndo.gPbo.sMnOa 


1 


111.51 


0.411 


-0.125 


125.8 


6 xlO -3 






±0.02 


±0.005 


±0.003 







The detailed results of fitting are given in Table 3. We have extracted the initial value of 7 from the y-axis intercept 
in the plot of C/T vs T 2 (inset of figure 6). To determine the magnon contribution to the specific heat, we extracted 
the spin stiffness constants D from the magnetization data [27| using Bloch T 3 / 2 law (see table 3). Subsequently, the 
corresponding values of P3/2 are calculated from D using the relation (3 3/2 = 0.113 Ra 3 (k B /D f/ 2 . Here, R = 8.314 
J K" 1 mole" 1 is the universal gas constant and a is the lattice parameter of the elementary perovskite cell [27[. Then 
the remaining coefficients B3 and B§ have been determined by fitting. The coefficient B3 corresponds to the Debye 
contribution to the specific heat at low temperature which can be expressed as, 

C D ebyc = (12/5) rRir 4 ( £-) (3) 
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where r is number of atoms in the unit cell, i.e. r = 5, R the universal gas constant, and 6jy is the Debye 
temperature [3^| . We have calculated the values of 6 D for x = 0.15, 0.3, 0.4 and 0.5 and they are 332 K, 292 K, 
291 K and 285 K, respectively. Generally, the electronic specific heat term is not expected to be present for samples 
with x = 0.15, 0.3 which are insulators at low temperature. However, it is observed that the magnitude of 7 is 
unusually large for all these samples. The enhanced values of 7 are most probably not related to the conduction 
electrons. The corresponding specific energy contribution may be due to magnetic effects related to the Nd ions and 
the Mn-sublattice [27]]. In particular, the magnitude of 7 is presumably influenced by the tail of the Schottky-like 
anomaly due to the presence of Nd ions. 

As can be seen from Figure 5, the specific heat below 15 K has a strong Schottky-like anomaly for all the samples. 
This effect is due to the Zeeman-like splitting of crystal- field ground state multiplets in the Nd 3+ ions, which was 
already described for Ndo.67Sr 33 Mn03 by Gordon et al. [35|. A similar phenomenon was reported for Pr .gSr .2MnO3 
[36] ]. Nd 3+ ions have a ten- fold degeneracy for the ground-state J multiplet 4 Ig/2 which is split by the crystal- field 
into five Kramers Doublets [37j • An effective molecular field H m f is assumed to be present at Nd sites. Although 
this field splits each of the five crystal-field doublets, at low temperature, only the ground state doublet needs to be 
considered 35]. Assuming that the effective moment of Nd 3+ ions in the ground state to be fj,Nd and that the splitting 
of the doublet is A s — 2fiNd H m f, a contribution from a two- level Schottky function should fit the excess specific heat 
at low temperature (T <15 K). We have used the Schottky function for a two-level system as (34l. I36T] 

C Sch (T,H) =n Sch N a k B 

where ns c h is the coefficient of the contribution from the Schottky effect, and N a is the Avogadro number. The fitting 
is carried out by adding this Schottky term to Equation (3), while keeping the lattice contribution (parameters B3, 
B 5 ) and the spin-wave contribution (^3/2) fixed. The result is shown in figure 7. The values of Schottky gaps A s 
and Schottky co-efficients ns c h for x = 0.15, 0.3, 0.4, 0.5 obtained by fitting are 1.14, 0.95, 0.88, 0.68 meV and 0.5, 
0.54, 0.56, 0.34 respectively. It has been noticed that this anomaly can be better fitted by a modification of the linear 
contribution 7T in the temperature range 2 to 15 K . For example, if we relax the value of 7 during fitting, the 
modified 7 values for x = 0.15, 0.3, 0.4 and 0.5 will become 44, 37, 29 and 61 mJ mole -1 K~ 2 respectively. Thus, 
the fact that the enhanced linear contribution 7T to the specific heat (as shown in table 2) is affected by magnetic 
contributions is reflected by these reduced 7 values [27| • We assume this broad linear specific heat contribution to 
originate from possible ordering of Nd moments. We observe that the magnitude of the Schottky gap A s increases 
with the concentration of Nd 3+ (figure 8) ions. Equivalently, the molecular field experienced by Nd ion is the strongest 
at the lowest x and decreases as x increases. However, the expected full contribution of the split ground-state Kramers 
doublet from the Nd 3+ -ions is not found in this fit, which would require that ns c h — 1 — x. Also the fit in the low 
temperature range is not overall satisfactory as seen in figure 7. 



(l + exp(^))2 



(4) 



IV. DISCUSSION 



We have seen that the gap energy A in the ferromagnetic insulating phase at low temperature varies with x 
and vanishes at a critical value x c w 0.33. Although the origin of the gap is not clearly understood yet, we can 
find some explanation in the light of the theory by Pai et al. (23J. The theory considers three important on-site 
interactions in manganites, namely the Jahn- Teller (JT) effect, Hund's rule coupling (Jh), and Coulomb repulsion U. 
The theory is based on the new idea of co-existing localized JT polaronic (I) and broad band (b) itinerant electrons 
from e g states. The effective bandwidth 2W (W — x 1 ' 2 Dq) of the b band decreases significantly as x decreases for 
any sizable U. Consequently, the bottom of the b band shifts above the Fermi level for small x. At T — all e g 
electrons become localized as I polarons. Mobile b states are occupied only by thermal excitations across the gap. 
The system is still ferromagnetic because of the Hund's rule coupling Jh that remains operative also for I polarons. 
This explains the insulating ferromagnetic behaviour at low doping with a thermally activated electronic transport 
[38| . Furthermore, with the increasing x, W increases and beyond a critical concentration x c , when the band width 
equals the JT-distortion energy, W(x) — Ejt, the low temperature state becomes a ferromagnetic metal. 

We have found the values for Do and Ejt of Ndi-^Pb^MnOa from our analysis. These are of acceptable order 
of magnitude, but appreciably smaller than the values estimated by Pai et al. . The smaller value for Do may be 
related to the decrease of half-bandwidth due to smaller cationic radius in Ndi-^Pb^MnOa which is comparatively 
smaller than in the wide-band systems like Lai-^Sr^MnOs [23J. However, the measured activation energy A(x) for 
the composition Ndo.ssPbo.isMnOs does not fit well with the dependence expected by the theory of Pai etal. [23| . 
This is possibly related to appreciable antiferromagnetic couplings giving rise to stronger spin-disorder scattering at 
this low doping close to the insulating antiferromagnetic phase. This could cause an increased apparent activation 
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energy for the crystal with x = 0.15. It is noticed that the metallic phase sets in above x = 0.3 and for compositions 
beyond this level we have found a clear metal-insulator transition at high temperatures. 

The empirical relationship between the coefficient A for the electronic contribution to resistivity and the coefficient 
7 of the electronic specific heat has been found by Kadowaki and Woods with A/y 2 ss 1 x 10~ 5 /i cm (mole K 
2 /mJ) 2 [14(. For the metallic composition x = 0.42 we can use the determined values for A and 7 to estimate this 
ratio. We find a Kadowaki- Woods ratio around 850 x 10~ 5 [i ft cm (mole K 2 /mJ) 2 . This means that the properties 
of the metal-like manganites are far from those expected for a normal metal. Since the Kadowaki- Woods empirical 
relationship is only valid for pure metals, it is not strictly applicable to the present system. However, considering 
the fact that samples with x > 0.4 are metallic at low temperature we can assume that there should be a positive 
contribution from e-e scattering. In fact, we have found that there is still a positive contribution from e-e interaction 
with conventional exponent oc T 2 in x — 0.42. But, this contribution is not prominent in x — 0.5 at HoH =0, 
possibly due to the resistivity minimum found for this composition. We have attributed this resistivity minimum to 
weak localization in a disordered correlated system (2lj . Resistivity minima at low temperature are also a typical 
characteristic of the Kondo effect [3!|. However, the resistivity minima due to a Kondo effect should disappear in 
relatively weak applied fields, which is not observed in the present case. Moreover, the manganite system at x = 0.5 
behaves as a homogeneous ferromagnetic metal and not like a metallic alloy with dilute magnetic impurities. Hence, 
the resistivity minimum observed here for the x = 0.5 single crystal should not be related to a Kondo effect. 

On the other hand, the enhanced weak localization contribution p t term in the presence of magnetic field is an 
expected effect. According to Ref. [2l|, the e-e interaction in the presence of strong disorder is special and the magnetic 
field enhances this interaction which results in an increased resistivity. From the resistivity in the a;=0.5 sample, we 
also find a reduction in the p p term by a magnetic field, which is consistent with a reduced two-magnon scattering. 
Here, one should remember that it has not been possible to distinguish the T 4 5 term due to two-magnon scattering 
from the T 5 electron-phonon interaction term. 

The weak localization and T 1 / 2 e-e interaction should be attributed to the static disorder in these manganites caused 
by the mixed A-site substitution of tri-valent rare-earth Nd 3+ by divalent ions (Pb 2+ ), which gives rise to some static 
electronic effects. The magnitude of disorder can be quantified by the A-site variance a 2 =< r\ > — < ta > 2 
where ta is average cationic radius [4J. The a 2 for x — 0.3 is 0.0101 [27j. Taking this value into consideration the 
disorder for x = 0.42 and 0.5 is around 0.0141 and 0.0168 respectively. The difference between the fitting of resistivity 
data for x = 0.42 and 0.5 (see table 1 and 2) is clearly observed. Hence, the disorder may not be large enough to 
cause a prominent effect of weak localization in x = 0.42, whereas its influence is strong in the sample with x — 
0.5. In addition, it should be mentioned for the specific Ndi-^Pb^MnOa system, that the half-doped system does 
not show indications of a charge-ordered/orbital-ordered and antiferromagnetic behavior. This is consistent with 
the proposed phase diagram of Ndi-^Pb^MnOs, where no indication of a charge-ordered/orbital-ordered phase was 
observed around 50:50 composition^. 

The analysis of the specific heat data points to some particular effects in the magnetic system across the series of 
crystals. It is seen in table 3 that the value of D increases with x, because the average cationic radius (< ta >) 
increases with x too. Since D is directly proportional to the exchange integral, the higher D value implies larger 
exchange coupling resulting in higher Tq [271 ] . The fact that the x = 0.5 sample has the highest Tq supports the 
above explanations. The Debye temperature 9d decreases with the increase of x here (see table 3). However, the 
Debye temperatures are smaller than usual values reported for manganites [41 ], This may be due to the fact that 
the analysis has not been carried out in the constant 9d region and the value is affected by systematic drifts (27l. [34|. 
The electronic contribution 7T to the specific heat is unusually large and present even in the insulating phase at low 
temperature. We believe that this broad specific heat contribution is related to a magnetic ordering of Nd moments, 
and/or to frustrated glassy ordering in the coupled Nd-Mn magnetic system owing to the dilution of the rare-earth 
A-site. Such an effect would be poorly emulated by an enhanced value of 7 over the temperature range up to about 
30 K. Hence, an exact evaluation of Kadowaki- Woods ratio is difficult for two specific reasons: (i) Even where the 
AT 2 term is clearly present as in x =0.42, the 7 from specific heat is overestimated by the magnetic contributions 
to the low temperature specific heat. Therefore, there is no clear way to estimate the Kadowaki- Woods ratio even in 
that case. In this respect, the true electronic 7 may be smaller, thus the calculated K-W ratio is a lower estimate, 
(ii) As our study has suggested for x=0.5 the AT 2 term is not discernible in some cases. Thus, the material does not 
behave like a metallic system at all, where such a term in resistivity would result from e-e scattering. 

The specific heat measurements at low temperature show a strong influence of a Schottky-like anomaly. The 
peaks due to Schottky-like anomaly are significantly broadened. This may the effect of a distribution of ground- 
state splittings on the Nd due to an inhomogeneous molecular field, e.g. related to a Nd-Nd magnetic exchange 
and the dilution on the A-site. In particular, there may be a greater number of Nd-sites with very low molecular 
field and correspondingly low Schottky-like contribution to specific heat. This can explain why only a part of the 
expected contribution of the two-level systems to the specific heat is found in our fitting. However, efforts to fit the 
Schottky-anomaly with more than one doublet energy to mimick a distribution of molecular fields were inconclusive 
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owing to the great number of necessary parameters in such fittings. The molecular field leading to the Zeeman-split 
levels of the ground state doublet in Nd 3+ includes also contributions from the Nd-Mn exchange couplings [4^|. A 
direct coupling with the Mn-0 subsystem may require some particular form of superexchange like 4f-3d coupling. If a 
Nd-Mn exchange mechanisms is to be the most important contribution, then an inhomogeneity of this molecular field 
would require an inhomogeneous ordering of the Mn-sublattice. On the other hand, the molecular field as measured 
by the magnitude of the effective single Schottky gap changes with Nd ion concentration. It is not obvious, why and 
how the interaction between Nd and Mn-0 subsystem could yield a decreasing molecular field with increasing doping 
x. In particular, the internal Mn-0 exchange becomes stronger with increasing x as seen from the increasing Curie- 
temperature up to x = 0.5. This indicates that the change of the effective molecular field H m f is influenced by the Nd 
content, which may be due to an exchange interaction between the Nd ions, e.g., by a long-range indirect exchange of 
the RKKY-type. In the case of Nd-Sr system it has been noticed that Nd-Nd interaction is small compared to Nd-Mn 
interaction [35|. Therefore, the detailed understanding of the magnetic effects involving the diluted Nd 3+ sublattice 
on the A-sites requires further theoretical analysis of possible magnetic exchange mechanisms in manganites. 



V. CONCLUSIONS 



The low temperature transport property and specific heat of Ndi-^Pb^MnOa have been analyzed for x = 0.15, 0.3 
, 0.4, and 0.5. In the considered temperature range, activated behaviour is found in resistivity data for samples, which 
are ferromagnetic insulators. This activation gap A disappears beyond x = 0.3, where the true zero-temperature MI 
transition takes place at an estimated critical concentration xp ps 0.33. A positive e-e scattering term AT 2 describes 
the leading resistivity effect in the low temperature transport data for the metallic sample x = 0.42. However, the 
validity of T 2 contribution to resistivity from e-e scattering used in earlier fitting formula is questioned here, as it is 
not discernible for the composition x = 0.5, which still shows overall metallic properties. A negative T 1 / 2 contribution 
from e-e scattering valid for disordered systems explains a shallow minimum in the resistivity for this x = 0.5 crystal. 
Low temperature specific heat data have been analyzed for these samples. Anomalously large values of the electronic 
specific heat coefficient 7 probably originate from magnetic effects. They could be related to a glassy magnetic 
ordering on the diluted Nd-sublattice. Indications of inhomogeneous magnetic properties are also found from a wide 
and subdued Schottky-like contribution to the specific heat from the Zeeman-split Kramers ground-state doublets in 
Nd 3+ -ions. 
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FIG. 1: The logarithmic plots of resistivity as a function of (1000/T) for x = 0.15, 0.3, 0.42 and 0.5. Inset: the logarithmic 
plot of resistivity as a function of (1000/T) for x = 0.5 in the low temperature range. 
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FIG. 2: The variation of the low temperature gap A in the FMI phase as a function of Pb concentration for Ndi-^Pb^MnOs 
for x = 0.15, 0.2, 0.25, and 0.3 (data of x — 0.2 and 0.25 are taken from Ref [Hf). The filled square points are experimental 
data. The dashed line is the result of a fit for the data x = 0.15-0.3. The thick solid line is the result of fit for x = 0.2-0.3, 
which is extrapolated to the critical concentration x c (thin line). 
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FIG. 3: (a) The resistivity data for x = 0.42 and results of fit according to Eq. (1) without the T 1 / 2 term. The inset shows the 
plot of resistivity data (after subtracting the residual resistivity) p vs T 1 ^ 2 in zero field. The points are experimental data and 
the solid lines are the fit. (b) The resistivity data for x = 0.5 at low temperature and the results of the fit according to Eq. (1) 
excluding the A T 2 term at and 5T. The points are experimental data and solid lines are fitting results. Inset shows the plots 
of resistivity data p vs T 1 / 2 for x = 0.5 at 0T and 5T. The data in both plots below 10 K show the presence of a negative T 1 ^ 2 
contribution due to e-e interaction for x= 0.5 sample. The straight lines are results of linear fit. 
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FIG. 4: The plot of A^ 1 dp/dT 2 vs T for x — 0.42, 0.5. This shows the deviation from the T 2 behaviour of the resistivity 
data at very low temperature for x = 0.5. The negative values of A -1 dp/dT 2 for x = 0.5 appear because there is a minimum 
and change of slope in the resistivity data. However, the data for the x= 0.42 sample follow the conventional T 2 behaviour. 
Inset(a) and (b) show p vs T 2 (after subtracting the residual resistivity) for x = 0.42 and 0.5 respectively. The straight lines 
are linear fits. 
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FIG. 6: Specific heat vs temperature data for the x = 0.4 sample. The solid line is the fit using Eq. (2) and open circles are 
experimental data. Inset shows C/T vs T 2 plot for x—OA 
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FIG. 7: Result of fits including the Schottky effect to the specific heat data of Ndi-^Pb^MnOs for x = 0.15, 0.3 and 0.4 in 
the low temperature range. The points are experimental and the solid lines are the fits. The data for x = 0.3 have been taken 
from Ref.[23] 




FIG. 8: The variation of the Schottky gap A s and the molecular field H m f with Pb concentration x for Ndi-zPb^MnOs. 



